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EFFECTS OF SOIL STRUCTURE ON* BURROW CHARACTERISTICS 
OF FIVE SMALL MAMMAL SPECIES 


John W. Laundré! and Timothy D. Reynolds? 


ABS [teAC 


Burrows of small mammals can impact a variety of soil processes including organic turnover, aeration, 


and anineralization rates. The strueture of burrows. depth, length. and complexity can influence the extent of the impact 
burrows have on soil proeesses. Soil properties. in turn, are thought to affeet burvew structure. To increase our under- 
standing of burrow-soil dynamics, we compared maximum depth, total volume, total length, volume:length ratio, and 
complesity of burrows of five small mammal species with bulk density and soil texture in multiple regression analyses. 
Burrows of Wyoming ground squirrels (Spermophilus elegans) were deeper, longer, and more complex as percentage of 
silt and clay increased and percentage of sand and bulk density decreased. Average maximum depth of montane vole 
Microtus montanus) burrows increased as soils became sandier. Length and volume of deer mice (Peromyscus manicu- 
latus. burrows inercased with increases in bulk density and percentage of clay. Volume. length. and complexity of kan- 
garoo rat) Dipodomys ordii) burrows were greater in soils with higher amounts of clay and silt. Townsend's ground 
squirrel Spermophilus townsendii) burrows did not appear to be affeeted by the soil properties measured. 


Key words: burrow structure. soil, small mammals, bulk density, soil texture, Idaho. 


The association between burrowing activi- 
tv of mammals and the subterranean environ- 
inent is an interdependent relationship that is 
receiving increasing attention because of its 
potential impact on plant community struc- 
ture (Andersen and MaeMahon 1955, Inouve 
etal. 19ST). Burrowing activity of mammals 
can impact a variety of soil processes inchud- 
ing organic matter turnover, inorganic distri- 
bution, aeration, and mineralization rates 
Abaturov 1972, Chew 1978. Zlotin and Kho- 
dashova 1980, Hole 1981, Reichman and 
Sinith 1989). The extent to which burrowing 
mammals may influence soil processes can 
depend on the strueture (complexity and 
dimensions) of burrow systems. Deeper bur- 
rows increase the depth, and more extensive, 
convoluted burrows increase the area of influ- 
ence compared with shallower and simpler 
systems. Thus, factors that influence burrow 
structure could ultimately determine the 
Hnpact a burrow has on soil processes in an 
area. 

Burrow structure varies considerably 
among species | Reynolds and Wakkinen 1987, 
Reichman and Smith 1989) and also differs 
within a species with length of oceupation of 
the burrow or age of occupant ‘Reichman and 
Smith 1959) Intraspecific differences in bur- 
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row dimensions are also hypothesized to be 
related to physical properties of soils, e.g., 
bulk density, texture, etc. (Anderson and 
Allred 1964, Revnolds and Wakkinen 1987). 
However, this hypothesis has not been ade- 
quately tested. To more clearly understand 
the dynamic relationship between burrowing 
mammals and the soil they live in, more infor- 
mation is needed on how soil characteristics 
nnpact burrow structure. Assuming soil prop- 
erties impact burrow structure, we predicted 
that burrow dimensions within a species 
should be quantitatively related to changes in 
measurable soil attributes. The purpose of this 
research was to test this prediction. 


STUDY AREA AND METHODS 


To test our prediction, data of Reynolds 
and Wakkinen (1987) were supplemented 
with data on burrows from several other soil 
types. All data came from within two study 
areas in southeastern Idaho: the Idaho 
National Engineering Laboratory (INEL), 65 
km NNW of Pocatello, Bannoek County, and 
public and private lands near Soda Springs. 
Caribou County. The INEL, a National Envi- 
ronmental Research Park on the upper Snake 
River plain, is classified as a cool sagebrush 


vorsity, Pocatello Tdaho $3209 
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(Artemisia spp.) desert dominated by sage- 
brush and grasses (Anderson and Holte 19S]. 
Sample sites near Soda Springs were in three 
mountain valleys: Wooley Valley, T33N. 
R43E, Sec. 25; Conda Valley. TSS, R4A2E, See. 
23; and Big Canyon, TOS, RA3E, Sec. 13. Vee- 
etation at these sites is also a sagebrush-grass 
mixture. 

Multiple regression analysis was used to 
test the prediction that burrow dimensions 
within a species are quantitatively related to 


changes in soil properties. Measurements of 


eaeh burrow charactcristic taken were indi- 
vidually used as the dependent variables: 
properties of the soils found at corresponding 
burrow sites were the independent variables. 
The null hypothesis tested for each analysis 
was that the multiple regression coefficient 
was equal to zero. 

For each species, burrows were excavated 
from at least four different subsites, repre- 
senting a variety of soil tvpes, within the 
study areas. Subsites were selected based on 
the presence of burrows in the area. Within 
the subsites, usually all the burrows were 
sampled. To quantify soil differences among 
burrow sites within a species, we measured 
bulk density (g/cm?) and soil texture (percent 
sand, silt, and clay) at approximately 10 cm 
below the soil surface at burrow locations 
with the core technique (Blake 1965) and the 
hydrometer technique (Day 1965), respective- 
ly. Soils in the study areas had relatively uni- 
form profiles with little vertical development, 
so near-surface measurements were consid- 
ercd adequate to classify the soil profile. [t is 
recognized that other soil properties exist and 
might influence burrowing behavior. Howev- 
er, estimates of density and texture are easily 
obtained, and other properties such as 
drainage, structure, and consistency are cor- 
related with these two measurements | Foth 
1978). 

The measurements we took to quantify the 
soils are highly correlated and would not be 
appropriate for use in a regression analysis. To 
classify soil types based on uncorrelated 
indices, we used the data on bulk density and 
soil texture to generate z-score standardized 
principal component scores (Manly 1990) for 
each soil sample taken at burrow sites. We 
then used the first two principal components 
(soil components) as the uncorrelated inde- 
pendent variables in our regression analyses. 
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All principal component analyses POA were 
conducted with the Biostat’ 3} Sigma Soft 
1130 Shalanwood Land Placentia. CA 92670 
statistical package. 

Burrows of five snd] maminal species 
were included in this study. Townsend's 
ground squirrel |Spermophilus townsendis 
Wrvoming vround squirrel |S. elegans). deer 
Ord s kan- 
Wand montane 
vole (Microtus montanus). Before we exam- 
ined a burrow system, we determined the 
species occupying the burrow by imuking 
visual observations, Snap-trappmg. or e\amin- 


mouse | Peromyscus maniculatuss. 
garoo rat) Dipadomys ordii 


ing hair and feces near burrow entrances. Lit- 
tle data were available on the age of burrow 
occupants. Of the animals taken with snap- 
traps. most were adults. We injected the bur- 
rows with polyurethane foam | Felthauser and 
MieInroy 1953), excavated the surrounding 
soil, and measured the burrow system. 

hive measurements were recorded: mani- 
mum depth. total length, vohime, volume. 
length ratio, and complexity. Maximum bur- 
row depth to the bottom of the burrow was 
determined by measuring depth in situ at [0- 
em intervals along the total length of the 
exposed foam casts. Total burrow length was 
the sum of the main and all side tunnels. VMter 
depth and length were recorded. casts were 
removed and burrow volumes estimated by 
water displacement |Revnolds and Wakkinen 
19S7, Revnolds and Laundré 1988). Volnme. 
length ratios were calculated and used as rela- 
tive indices of burrow diameter to distineuish 
among different length-diameter ratios. c.g. 
long-narrow vs. short-wide burrosss within a 
species. The complesity of burrows was caleu- 
lated as the length of a dine connecting the 
two most distant points of the burrow divided 
by the burrow’s total length. This index is 1.0 
for linear burrows and progressively Tess with 
increased burrow complexity. This techmque 
was found superior to other “indices of tnear- 
ity ‘Reichman et al. 1952. Cameron etal 199 
that approached infinity for the short. rela 
tively linear burrow systems found in thas 
study. Additionally, the five burross character 
istics measured were sed to Generate z-score 
standardized principal component scores for 
the burrows of cach species. The first priicipal 
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\ \Ieans ) f micasurements + SD) recorded for burrows examined in this study. The lower portion of the 
tele sontains eizenvector coefficients for maximum depth cm), volume (1), length (m), volume per length (vol/len), and 
complexity for the first /Z1-b principal component from the principal component analysis of burrow characteristics. 
The percent variance var explained by the principal component is also given. 
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“Deennouse Peromyscus maniculatus 
 Kangarov tat) Dipodomys ordit 
Montin 
‘Townsend's ground squirrel Spermophilus townsendit 


vole Mierotus montanus 


Wyonung ground squirrel Spermophilus elegans 


component score (burrow component, 
71-b) of a burrow was used to characterize a 
burrow system as a whole. Burrow compo- 
nents within a species were then used as the 
dependent variable in a multiple regression 
analysis with the first and second soil compo- 
nents. This analysis would help determine 
whether the burrow system as a whole was 
influenced by soil properties. 

Ground squirrels construct “shallow” and 
“deep” burrows (Revnolds and Wakkinen 
1987). To test for relationships between shal- 
low or deep systems and soil types, we reana- 
Ivzed separately the maximum depth of the 
two system types for the two ground squirrel 
species. Shallow burrows for Townsend's 
ground squirrels were classified as being less 
than 60 cm deep (Reynolds and Wakkinen 
1987). Based on an analysis of frequency dis- 
tributions of maximum depth found in this 
study (unpublished datas. Wyoming ground 
squirrcl burrows less than 90 cm deep were 
classified as shallow systems. After the bur- 
rows were reclassified as shallow or deep, we 
found we had sufficient sample sizes to con- 
duet only our regression analyses for shallow 
burrow Ys 

Statistical comparisons were calculated 
with Statistics with Finesse® (J. T. Bolding, 
Box 339, Favetteville. AR 72702), Biostatis- 
tics® Lin cotts Directory. 40 Glen Drive, 
Mill Valley, CA 94941. or Biostat® (Sigma 


Soft, 1430 Shalanwood Lane, Placentia, CA 
92670) computer packages. All percent data 
were arcsine transformed before statistical] 
calculations were made. The significance 
level was P = .05, and all reported means are 
+ standard deviation. 


RESULTS 


One hundred forty-nine burrow systems 
were excavated. Means and ranges of mea- 
surements recorded are presented in Table 1. 
Included in Table 1 are eigenvector coeffi- 
cients for the first principal burrow compo- 
nents (Z1-b) for the five species. Except for 
montane voles, the first principal components 
explained 52-59% of the variance of the bur- 
rows sampled. 

Bulk density and texture of soils at burrow 
sites varied among the five species and were 
quite variable within a species (Table 2). The 
first (Z1-s) and second (Z2-s) principal compo- 
nents of soils from burrow sites explained 
over 90% of the variance in samples for all 
five species (Table 2). Eigenvector coefficients 
for the equations describing the relationship 
among the soil properties for the first and sec- 
ond principal soil components are also pre- 
sented in Table 1. For the first principal soil 
component equation within a species, bulk 
density and percent sand had the same sign 
(+ or -), which was opposite the signs for per- 
cent silt and clay. 
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TABLe 2. Means ‘+ SD) aud ranges of bulk density em? and texture §¢ 
burrow sites. Testure values are untransformed. The lower portion of the table conlans eng 
bulk density (BD) and percentages of sand. st. and clay for the first (Z]-s . 
from the principal component analyses PCA) of burrow site soils 


cipal component is also given. 


Species Bulk density: “ sand 
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“Deermouse | Peromyscus maniculatus 

bKangaroo rat/ Dipodomys ordin 

Montane vole) Microtus montanus 

Townsend's ground squirrel: Spermophilus townsendir 
Wyoming ground squirrel] Spermophilus elegans 


Of the 28 multiple regression analyses of 
burrow characteristics with the first and see- 
ond principal soil components, five were sig- 
nificant (Table 3). Four of the five significant 
multiple regressions were for characteristics 
of Wyoming ground squirrel burrows. Coeffi- 
cients of determination for these four signifi- 
cant multiple regressions ranged from 17 for 
maximum depth of shallow burrows to 27% 
for maximum depth of all burrows combined 
(total depth). For these five significant imulti- 
ple regressions, the only partial regression 
coefficient for the first principal soil compo- 
nent that was not significant was for total max- 
imum depth (Table 3). In addition to signifi- 
cant multiple regressions, partial regression 
coefficients of the first prineipal soil compo- 
nents were significant for volume per length 
and complexity. Partial regression coefficients 
of the second principal soil components were 
significant only for total maximum depth, 
shallow maximum depth, and length of bur- 
rows (Table 3). 


For Wyoming ground squirrels. the multi- 
ple regression coefficient for the first princi- 
pal burrow component Z)-b) was alse signifi- 
cant (Table 3). Of the two independent vari- 
ables, only the first principal soil component 
Z1-s) had a significant partial regression coel- 
ficient. There was an inverse reletorship 
between Zl-b and Z]-s. and the regression 
equation explained 16° of the variability: seen 
in Zl-btFig. Hh. 

The fifth significant multiple regression 
was for depth of taontane vole burraws ‘The 
coefficient of determination for this riultiple 
Table 3). Partial recres- 
sion coefficients far both independent v.a- 
Zles and Z2-s) were sienifieant and 


regression was 16 


ables 
negative. The regression equation for depth 
on Zl-s explained most | 106) of the variabili- 
ty seen in depth of microtus burrows Fig. 2 

The partial regression coeficient of the first 
principal soil component was also significant 
for the volume-length ratios of microtus bur- 
rows. The regression coefficient was positive, 
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Taper 3. Results of multiple linear regression analyses between burrow characteristics, including the first principal 
burrow component Zi1)], and the first two principal components for burrow site soils for the five species examined. 
Only regressions that were statistically significant are presented. Coeflicient of determination (r?), results of signifi- 
cance tests F), and probability (P) are given for cach multiple regression analysis. Significance values of partial regres- 
sion coefficients for the principal soil components are also presented. Signs in front of f values indicate direction of the 
regression relationship. 
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ume and length of burrows and for the first 
principal burrow component (Table 3). Coeffi- 
cients of determination were 0.23 for volume 
(Fig. 4a) and 0.22 for length (Fig. 4b). The 
regression equation for Z1l-b on Z2-s 
explained 38% of the variability in the data 


and the regression equation explained 11% of 
the variability in the data (Fig 2b). 

None of the multiple regression eoeffi- 
cients for burrow characteristics were signifi- 
cant for kangaroo rats, deer mice, or 
Townsend's ground squirrels. However, par- 


tial regression cocllicients for the first prinei- 
pal soil components were significant and posi- 
tive for volume and length of kangaroo rat 
burrows (Table 3). Regression equations 
explained 27% of the variability in volume 


‘Hig. 3a and 36% of the variability in length of 


kangaroo rat burrows (Fig. 3b). Partial regres- 
sion cocticicents for the first prineipal soil 
component: Z]}-s) and the first prineipal bur- 
row component (Z1-b) for this species were 
also significant but negative. The regression 


equation for Zi-b on ZI-s explained 25% of 
the variab n the data (Fig. 3e). 
or deei artial regression coeffi- 


eients for the s 


‘ipal soil component 
(72-s) were significan 


yositive for vol- 


(Fig. te). 

For Townsend's ground squirrels, the only 
significant partial regression coeffietent was 
for burrow volume on the second principal 
soil component (Table 3); the regression equa- 
tion explained 28% of the variability (Fig. 1). 


DISCUSSION 


Our results indicate that differences in soil 
properties have different effects on the 
species of small mammals studied. Town- 
send’s ground squirrels seemed the least 
affected by soil differences. For this species, 
only volume was significantly related to the 
second principal soil component. The second 
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Wyoming Ground Squirrel 


Y = 0.03 ~ 0.45 X 
r 


First PC (Z1) scores for burrows 





First PC (Z1) scores for soils 


Townsend's Ground Squirrel 


= 12.14 + 9.52 X 
r = 028 


Burrow volume (1) 





Second PC (Z2) score for soils 


Fig. 1. Regressions, with 95% confidence intervals. of 
the first principal barrow component scores on the first 
principal soil component scores for Wyoming ground 
squirre]) burrows, and of burrow volime on the second 
principal soi] component scores for Townsend's ground 
squirrel burrows. All principal component scores are z- 
score standardized. 


principal soil component was influenced pasi- 
tively by bulk density and percent silt and 
negatively by percent sand (Table 2). Thus. 
Townsend's ground squirrels constructed larg- 
er burrows in the firmer loamy soils. Reynolds 
and Wakkinen (1987) found significant regres- 
sion relationships for this species among, soil 
texture components and burrow depth. 
length, and vohime. However, they used indi- 
vidual soi] separates (percent sand, silt, and 
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Fig. 2. Regressions. with 95% conlidence mternals oo 


burrow depth a and volume per length vol lem rats 


b oon the first principal soi component scores for myo 
tane vole burrows. AT) prinerpal component scores grt 
score standardized. 


clay) in their analysis. As percent sand. silt. 
and clay are highly correlated. them ase could 
bias the results of a regression analysis and 
account for the different findings between the 
two studies. 

for kangaroo rats, burrow length and vol- 
ume seemed directh influenced by the per- 
centage of siltand clay in the soil. There was a 
significant negative regression of the first 
principal burrow component .Z1-by on the 
first principal soil component | Z1-s)) Table 3 
The magnitude of the first principal sail cam- 
ponent is directly related to percent silt and 
percent clay and inversely related to percent 
sand and bulk density. Table 2). The first prin- 
cipal burrow component is inversely related 
to all burrow measnrements except compleri- 
ty (Lable 1. Hawever a deerease in the com- 
plexity index indicates more complex hur- 
rows. Of the measurements, burrow depth 
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Fig. 3. Regressions, with 93° confidence intervals. of 
volume ‘a, length (b). and first principal burrow compo- 
nent scores (c) on the first principal soil component 
scores for kangaroo rat burrows. All principal component 
scores are z-score standardized. 


contributes little. Consequently, as the per- 
centages of silt and clay increase (high Z1-s). 
burrows do not become much deeper but 
thes do become longer and more complex 

“ complexity index). resulting in greater 


volyincs of soil being removed (low Z1-b). 
This interpretation is supported by the signifi- 
Cant and positive partial regression coeffi- 
Cie und for burrow volume and length 
with | rincipal soi] component (Table 
3). R ind Wakkinen (1957) did not 
find a : miship between bur- 
row cha iG soil texture for this 
species. Again I .in data analysis 


likely accounted for thei diff findings. 
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Second PC (Z2) scores for soils 
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Second PC (Z2) scores for soils 


Y = 1.08X - 0.23 
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Second PC (Z2) scores for soils 


Fig. 4. Regressions, with 95% confidence intervals, of 
volume (a), length (b), and first principal burrow compo- 
nent scores {c) on the second principal soil component 
scores for deer mice burrows. All principal component 
scores are z-score standardized. 


Several deer mice burrow characteristics, 
ineluding the first principal burrow compo- 
nent, seemed to be influenced by soil proper- 
ties as expressed by the second rather than 
the first principal soil component. For deer 
mice, the magnitude of the second principal 
soil component (Z2-s) increases as percent 
clay and bulk density increase or percent silt 
decreases; percent sand does not appear to 
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influence this second soil component |Table 
2). Relative to the first principal burrow com- 
ponent (Z1-b), all burrow measurements 
except complexity increased with higher com- 
ponent scores (Table 1). However, volume 
length measurements did not contribute sub- 
stantially to the score. Consequently. deer 
mice constructed deeper, longer, and more 
complex burrows (high Z1]-b) in soils with 
higher bulk densities and percent clay high 
Z2-S). 

Burrowing behaviors of Wyoming ground 
squirrels seemed the most affected by differ- 
ences in soil properties. Wyoming ground 
squirrel burrows had the most significant 
regression coefficients between characteris- 
tics and principal soi] components. The first 
principal soil component (ZI-s) increased 
with increasing silt and clay and decreased 
with increasing bulk density and sand (Table 
2). The first principal burrow component Z1- 
b) was negatively related to all burrow prop- 
erties except complexity (Table 1). As percent 
silt and clay increase (high Z1-s). all burrow 
characteristics increased (low Z1-)). 

For montane voles. only depth and vol- 
ume/length ratios of their burrows secmed 
influenced by soil properties. The first princi- 
pal soil component was positively influenced 
by percent silt and clay and negatively related 
to percent sand and bulk density |Table 2). 
The relationship between depth and the first 
principal soil component was negative. Con- 
versely, the volume/length ratio was directly 
related to the first principal soil component. 
Consequently, voles constructed deeper but 
narrower burrows in sandier but firmer soils. 

In summary, five measurements supported 
the prediction that burrow characteristics are 
affected by soil properties. Most of these five 
characteristics were of Wyoming ground 
squirrel burrows. and we conclude that the 
burrows of this species are influenced in a 
predictable manner by the soil properties 
measured. Ten of the remaining burrow mea- 
surements were also significantly influenced 
by soil properties as described by individual 
principal soil components. We conclude that 
the remaining burrow measurements were 


not affected by soil attributes. Variability of 


these burrow properties is likely influenced 
by other, yet to be determined, factors such as 
length of occupancy or age and sex of occu- 
pant. To investigate the influence these factors 
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may have on burrow structure and further 
delineate the WMpact of som propertres we 
suggest that controlled experiments be con- 
ducted. 

Our results suggest that soil characteristics 
of an area can affect various burrow dimien- 
sions in a predictive manner for the five sna 
mammal species we studied. Soil effects on 
burrow structure could. in turn. mfluence the 
soil processes of that area. A difference im 
maximum depth of burrows changes the loca- 
tion of the reservoir of nutrients for reevcling. 
increases the depth of soil aeration. and. espe- 
cially in arid and semiarid arcas. alters shal- 
low subsurface water recharge patterns. Dil- 
ferences in volume. length. volume.length 
ratio. and complexity probably have a greater 
impact on the magnitude rather than the 
direction of the influence of a burrow an soil 
properties. For burrows of the same species 
with the same maximum depth. longer, larger. 
or more coinples burrows within a given area 
would result in more surface soil deposition 
from. more aeration of. and more water infil- 
tration to a given profile depth than shorter. 
smaller. or simpler burrows. 
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